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Abstract. Any infinite sequence of substitutions with the same matrix 
of the Pisot type defines a symbolic dynamical system which is minimal. 
We prove that, to any such sequence, we can associate a compact set 
(Rauzy fractal) by projection of the stepped line associated with an 
element of the symbolic system on the contracting space of the matrix. 
We show that this Rauzy fractal depends continuously on the sequence 
of substitutions, and investigate some of its properties; in some cases, 
this construction gives a geometric model for the symbolic dynamical 
system. 



1. Introduction 

It is well-known that any unimodular irreducible Pisot substitution defines 
a dynamical system that is a finite extension of a toral translation, and it is 
conjectured that such a dynamical system is in fact measurably equivalent to 
a toral translation, see [7], [H]- A geometric model of the symbolic system 
can be obtained by projecting the discrete line associated to a fixed point of 
the substitution along its asymptotic direction, as we explain more precisely 
in the next section. 

In some cases, this can be extended to systems generated by an infinite se- 
quence of substitutions belonging to a finite set S (so-called 5-adic systems), 
the best example being that of the sturmian sequences, which are almost 1-1 
extensions of irrational circle rotations, the adic expansion being given by 
the additive continuous fraction expansion of the angle. 

It would be very interesting to be able to generalize this property; however, 
in the general case, we cannot expect it to hold for all sequences: already, 
one can check that the closure of the set of all sturmian sequences contains 
periodic sequences, with bounded complexity and finite 5-adic expansion; it 
is a degenerate case, where the symbolic model is finite. The existence of non 
balanced episturmian sequences (see |5]) is another obstruction; in that case, 
the classical construction of the Rauzy fractal by projection cannot work, 
since the projected set is well defined, but not bounded. The existence of 
minimal, but not uniquely ergodic S-adic systems associated with interval 
exchange maps is a third one : for the corresponding symbolic system, the 
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frequency is not denned, hence the discrete line (see below) associated with 
a symbolic sequence has no asymptotic direction, so the projection is not 
even defined. 

In this paper, we solve the problem in a restricted case: we consider a 
matrix A with positive integer coefficients of the Pisot type (that is, all 
eigenvalues are nonzero, and all eigenvalues except one are strictly smaller 
than 1 in modulus), and a finite set S = {<7i, . . . , <7fc} of substitutions with 
the same matrix A. We first prove that any infinite sequence of elements of 
S defines a minimal symbolic system, by defining a generalized fixed point 
(the limit point). We then prove, using a generalized prefix-suffix expansion, 
that this limit point stays within bounded distance of the expanding line of 
the matrix, and that we can associate a compact set by projection of this 
limit point on the contracting space. We also prove, by using a generaliza- 
tion of the classical IFS theorem of Hutchinson [9j, that this set depends 
continuously on the sequence of substitutions, and investigate some of its 
properties. 

In section 2, we fix the notations, and show that a primitive sequence 
of substitutions on an alphabet A defines a minimal dynamical system in 
„4 N . In section 3, we prove that if all the substitutions have the same inci- 
dence matrix of the Pisot type, we can generalize to that case the projection 
construction of the Rauzy fractal. In section 4, we generalize the classical 
iterated function system to the case of an infinite sequence of contractions, 
and apply it to our case to prove that the generalized Rauzy fractal depends 
continuously on the sequence of substitutions. In section 5, we give a few 
examples, and make some remarks. 

2. Substitutions and adic systems 

2.1. General setting. Let A := {1, d} be a finite set of cardinality d > 1, 
called the alphabet. We denote by A* the free monoid on A, (set of finite 
words on the alphabet A, with empty word denoted by s, endowed with the 
concatenation map). We denote by A N the set of infinite sequences on A, 
with the natural product topology. 

The length \w\ of a word w S A n with n £ N is defined as \w\ = n. 
For any letter a € A, we denote the number of occurrences of a in w = 
w\W2 ■ ■ -Wn-iWn by \w\ a = ${i\wi = a}. We will denote by I : A* h- > 
: w I—?- (|t«| a ) ag ^ 6 N d the natural homomorphism (abelianization map) 
obtained by abelianization of the free monoid. 

A substitution over the alphabet A is a nonerasing endomorphism of the 
free monoid ^4*. To any substitution a, one can associate its incidence matrix 
M, which is the d x d matrix obtained by abelianization, i.e. Mij = 
By construction, one has l(a(w)) = Ml(w) for any word w S ^4*. 

Definition 2.1. A substitution a is primitive if there exists an integer k 
such that, for each pair (a,b) G A 2 , \a k (a)\b > 0. 
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It is equivalent to suppose that the incidence matrix is primitive, that is, 
this matrix has a strictly positive power. We will always suppose that the 
substitution is primitive, this implies that for all letter j £ A the length of 
the successive iterations o~ k {j) tends to infinity. 

Remark 1. Since the incidence matrix of a primitive substitution is a prim- 
itive matrix, by the Perron-Frobenius theorem, it has a simple real positive 
dominant eigenvalue (3. 

2.2. Dynamical system defined by a primitive substitution. A sub- 
stitution a naturally extends to the set A^ of infinite sequences with value in 
A, by defining a(u\U2 ■ ■ •) = a(ui)a(u2) ■ • •■ We say that a sequence u 6 A^ 
is a periodic point of a if there exists some integer k such that a k (u) = u. 
One easily proves that any primitive substitution has periodic points, and 
that two periodic points with the same initial letters are equal, hence there 
are at most d periodic points. 

We denote by S the shift on A^, defined by S(u) = v, where v is the 
sequence such that, for all i £ N, V{ = tij+i. 

Definition 2.2. Let a be a primitive substitution, and let u be a periodic 
point of a. Let Vt a = {S n (u)\n E N} be the closure of the orbit of u by the 
shift. 

The dynamical system defined by the substitution a is (Q a ,S). 

Remark 2. The primitivity condition implies that any word occuring in a 
periodic word of a also appears in any other periodic word of a, hence all 
these periodic words have same closure for the shift, and the set fl a does 
not depend on the particular periodic word we have chosen, only on the 
substitution a. 

2.3. Pisot substitutions. We recall a classical definition of algebraic num- 
ber theory: 

Definition 2.3. A Pisot number is an algebraic integer /3 > 1 such that each 
Galois conjugate /3W of (3 satisfies \ /3® \< 1. 

By analogy, we define substitutions of Pisot type and irreducible substi- 
tutions of Pisot type: 

Definition 2.4. A substitution is of Pisot type if the dominant eigenvalue 
of its incidence matrix is a Pisot number. A substitution of Pisot type is 
(algebraically) irreducible if the characteristic polynomial of its incidence 
matrix is irreducible over the field of the rational numbers. 

Remark 3. It is equivalent to say that all the eigenvalues of the incidence 
matrix are not equal to 0, and all except one are strictly smaller than 1 in 
modulus. Any irreducible substitution of Pisot type is primitive (see [3]), 
but the converse is not true: there are primitive substitutions of Pisot type 
which are reducible. 
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Remark 4. The fixed point of a substitution of Pisot type can never be 
periodic. 

2.4. »S-adic systems. Dynamical system defined by substitutions form a 
small class; in particular, there is only a countable number of such system. 
For example, the only circle rotations whose natural coding is given by a sub- 
sitution are those whose rotation number is a quadratic integer of a special 
type (Galois integer). 

It is however possible to enlarge this class by replacing one iterated sub- 
stitution by an infinite sequence of substitutions chosen in a finite set S; this 
can be considered as a kind of continued fraction expansion, the substitution 
dynamical systems corresponding to periodic continued fraction expansions. 
We obtain in this way the class of S'-adic systems, but some care must be 
taken to generalize the primitivity condition. 

Definition 2.5. Let S be a finite set of substitutions on a fixed alphabet A. 
Let (cr n )neN be a sequence of substitutions in S. We say that the sequence is 
primitive if, for any n £ N, there is an integer k such that the substitution 
cr n o~ n +i . . . cr n+ k has a strictly positive incidence matrix. 

We can extend the notion of periodic points to such a sequence of substi- 
tutions. 

Definition 2.6. Let (o- n )n£N be a sequence of substitutions in S. We say 
that u £ A^ is a limit point of the sequence (c ra ) ne N if there exists a sequence 
of points of A N such that u = and u^> = a n (u^ n+1 ^). 

Proposition 2.1. Any primitive sequence of substitutions has a finite and 
nonzero number of limit points. 

Any finite word that occurs in a limit point of a primitive sequence of 
substitutions occurs in all the other limit points. 

Proof. With any substitution a, we associate the map / : A — > A which 
sends any letter a 6 A to the first letter of a{a). Given a sequence (a n ) of 
substitutions, we consider the sequence /i(/2 • • • (f n (A)) . . .)) of subsets of A. 
This is a decreasing sequence of nonempty subsets of A, so their intersection 
is nonempty. Let a be any element of this intersection; by construction, there 
exists a sequence (a n ) of elements of A such that ao = a and a n = f n (a n +i) 
for all n. 

Consider the sequence of words U n = a^a\ . . . a n -i(a n ). Since a n is the 
first letter of o~ n (a n +i), U n is by construction a prefix of ?7 n +i. The primi- 
tivity condition implies that the length of U n tends to infinity, since for any 
n there is a k such that a n a n+ \ . . . a n+ k sends any letter to a word of length 
at least d > 1. 

This sequence of finite words defines a unique infinite sequence u which 
admits these words for prefix. By shifting the sequences (o~ n ) and (a n ), we 
define in the same way points which, by construction, satisfy = 
a n (u^ n+l ">). Hence any primitive sequence of substitutions has a limit point. 
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Remark that any limit point arises in this way: if we have a sequence 
(«(")) of points of A m such that «( n ) = a n {u (n+l ^), and if we define a n to be 
the first letter of u^ n \ the words co°i • • • cr n-i( a n) are by definition prefixes 
of u^°\ whose lengths tend to infinity. 

Since, for any sequence (a n ) such that a n = f n (a n+ i), a n defines a p for all 
p < n, there can be at most d such sequences (because there are at most d 
prefixes of length n of such sequences for any n). Hence there are at most d 
limit points. 

Let u and v be two limit points, and let (a n ) and (b n ) be the corresponding 
sequences of letters. Let U be a finite word which occurs in u. Then there 
exists an n such that U occurs in the word ctoUi . . . o n -\(a n ). By primitivity, 
we can find k such that the letter a n occurs in the word a n . . . a n+ k{b n+ k). 
Hence U occurs in the word c7o°"i • • • &n+k(b n +k)i which is a prefix of v. This 
proves that a word which occurs in a limit word occurs in all the other limit 
words. □ 

Remark 5. The special case of a constant sequence corresponds to the fixed 
point of a substitution; in that case, the definition of primitivity of a sequence 
reduces to the definition of primitivity of a substitution. 

Definition 2.7. Let S be a finite set of substitutions, and let a = (cr n ) ne N 
be a primitive sequence of substitutions. The S-adic system generated by the 
sequence of substitutions a is the dynamical system (Q a , S), where is the 
closure of the orbit by the shift of a limit point of a. 

The previous proposition shows that the set Q, a does not depend on the 
choice of the limit point u, but only on the sequence a. 

Proposition 2.2. The S-adic system generated by a primitive sequence of 
substitutions is minimal. 

Proof. It is enough to prove that any word which occurs in a point of the 
system occurs with bounded gap, and it is enough to prove this for a limit 
point. Using the same notations as above, let U be the finite word under 
consideration. Then we can find n such that U occurs in <7o°i • • • &n—i( a n)> 
and by primitivity we can find k such that U occurs in ctq . . .cr n +fc(6) for 
any b in A. Let L be the maximum length of the words uq . . . a n+ k{b); since 
the limit point can be decomposed in these words, the word U occurs in the 
limit point with gaps bounded by 2L. □ 

Example 2.1. A very simple example of a family of S'-adic system is given 
by the sturmian sequences, or rotation sequences, which are almost one-to- 
one extension of rotations. Here, the alphabet A = {0, 1}, and the set S of 
substitutions is S = {ao, ai}, where <ro is defined by o"o(0) = 0, <ro(l) = 10 
and <7i is defined by 0i(O) = 01, ci(l) = 1. In that case, a sequence is 
primitive if and only if it is not eventually constant, and the sequence of 
substitutions which defines a given sturmian system is given by the additive 
continued fraction expansion of the angle of the associated rotation. 
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Example 2.2. A more tricky example is given by the natural symbolic dynam- 
ics of interval exchange maps. Here also we get a finite set of substitutions, 
associated with the Rauzy induction which generalizes the classical contin- 
ued fraction. This example shows that we cannot expect to generalize to the 
S-adic case all the properties of substitution dynamical systems; in particu- 
lar, there are S-adic systems defined by a primitive sequence of substitutions 
which are not uniquely ergodic, as shown by the existence of minimal but 
not uniquely ergodic interval exchange maps. 

3. Geometric models for substitution dynamical systems and 

5-adic systems 

Substitution dynamical systems are one of the simplest families of dynam- 
ical systems, and they are present in every framework where self-similarity 
occurs. It is natural to study their dynamical properties; as we have seen 
above, they are minimal, and a simple application of the Perron-Frobenius 
theorem shows that any letter, and any finite word, has a well defined fre- 
quency, which proves that the system is uniquely ergodic. 

Hence, a substitution dynamical system is a measured dynamical system 
in a natural way, and we can study its ergodic properties, look for eigenvalues, 
and look for the maximal equicontinuous factor. It is well-known that the 
Pisot condition plays a role here: the work of Rauzy and other people 
[7] has shown that any unimodular irreducible substitution of Pisot type on 
d letters is a finite extension of a translation of the torus T d ~ l , and it is 
conjectured that it is an almost one-to-one extension. 

We will briefly review a geometric construction for this maximal equicon- 
tinuous factor of a unimodular Pisot substitution a. We first need to intro- 
duce some algebraic formalism in order to embed the periodic point u in a 
linear subspace spanned by the algebric conjugates of the dominant eigen- 
value of the incidence matrix of a; the closure of the "projections" of the 
prefixes of u will form the so-called central tile or Rauzy fractal. 

3.1. Rauzy fractal associated to a unimodular Pisot substitution. 

Definition 3.1. A stepped line L = (x n ) in M. d is a sequence (finite or 
infinite) of points in ~R d such that x n +i — x n belongs to a finite set. 

A canonical stepped line is a stepped line such that xq = and for all 
n->0, x n+ \ — x n belongs to the canonical basis ofM. d . 

Using the abelianization map, to any finite or infinite word W, we can 
associate a canonical stepped line in lR. d as the sequence (l(P n )), where P n 
is the prefix of length n of W. 

The Perron-Frobenius Theorem implies that the canonical stepped line 
associated with a periodic point of a primitive substitution has an asymptotic 
direction. If this substitution is irreducible and of Pisot type, this canonical 
stepped line remains within bounded distance from the expanding line (given 
by the Perron-Frobenius right eigenvector of M a ). These two properties 
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(bounded distance and Pisot type) are equivalent for algebraically irreducible 
substitutions. 

Notation 3.1. Let a be an irreducible Pisot substitution. We denote by 
E s the stable space (or contracting space) and E u the unstable space (or 
expanding line) of M a . We denote by tt s the linear projection on E s , parallel 
to E u . 

Since the canonical stepped line stays within bounded distance of E u , its 
projection is a bounded set in the (d — l)-dimensional vector space E s . 

Definition 3.2. Let a an be irreducible Pisot substitution. The Rauzy frac- 
tal (or central tile) associated with a is the closure of the projection of the 
canonical stepped line associated to any periodic point of a in the contracting 
plane parallel to the expanding direction. We will denote it by lZ a . 

The Rauzy fractal can be decomposed in subfiles depending on the letter 
associated to the vertex of the stepped line that is projected. For any i £ A, 
we define : 1Z l a := {-ir s (l(uo . . . Uk-i), k 6 N, Uk = i}. We have obviously 

Remark: One can prove that the definition of 1Z a and 1Z l a for i G A does 
not depend on the choice of the periodic point u £ A. 

We will briefly recall some properties of the Rauzy fractal. If we denote 
by T the projection by tt s of the diagonal lattice of Z rf (that is, the set of 
(ni, . . . , rid) such that J2 n i = 0), the translate by T of the Rauzy fractal 
cover the stable plane. Hence the Rauzy fractal has positive measure. The 
projection from the orbit of the periodic point to the Rauzy fractal extends 
by continuity to all of £l a ; if we quotient by T, this projection gives a semi- 
conjugacy between (Q a ,S) and the translation by ir s (e) on lZ a /T, where e 
is any vector in the canonical base (all ir s (e) are equivalent modulo T). 

It is conjectured that the semi-conjugacy is a conjugacy, TZ a is a funda- 
mental domain for T and the 1Z l a are a partition of 1Z a , up to a set of measure 
0, and this can be proved in many particular cases. 

Finally, the subsets of the Rauzy fractal are solution of a set equation: 

Proposition 3.1. Let a be a unimodular irreducible substitution of Pisot 
type, with incidence matrix M. Let us denote by li = \cr(i)\ the length of 

cr(i), and o(i) = . . . wj^ = P^w^ where pff* is the prefix of length 

k — 1 of o~{i) and is the suffix of length li — k of cr(i). 

Then the subsets 1Z l a of the Rauzy fractal satisfy the following set equation: 

K= |J M.TZi + nMPh) 

j,k;wl=i 

3.2. Rauzy fractal associated to a sequence of substitutions. We 

want to generalize this construction to the S'-adic case. To do this, we 
can replace the periodic point by a limit point, but we need to have an 
asymptotic direction, and a bounded distance property. The most general 
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condition under which this is possible is very unclear, and we will restrict to 
a very special case. 

We consider a positive integral primitive matrix M with determinant ±1 
satisfying the Pisot condition: all the eigenvalues different from the Perron- 
Frobenius eigenvalue are strictly smaller than 1 in modulus. Let S be a set 
of substitutions, all of them having M as incidence matrix. 

It is then clear by construction that any sequence in S* N is a primitive 
sequence of substitutions. To such a sequence, one can associate a limit 
point, as we have shown before. This limit point has an asymptotic direction, 
which is the unstable direction of the matrix M. 

We need to fix some notation, a = (<J n )neN £ <5 N is a sequence of substi- 
tutions; as above, we write, for any substitution do £ S, k = \&o(i)\ (it only 
depends on M, not a ), and a (i) = . ..w^ = P™'*™?'* 'S° '\ where 
P^ ' 1 is the prefix of length k — 1 of ao(i). 

Lemma 3.1. Let u be the limit point of a sequence a € S N . The canonical 
stepped line of u stays within bounded distance of the expanding line of M. 

Proof. Let U be a prefix of the limit point. By construction, we can write 
U = (j ([/ (1) )P , where is some prefix of a limit point of the shifted 
sequence Sa, and Pq is a prefix of ao(i) for some i G A. We can iterate on 
U^, and a recurrence proves that there is some k such that we can write 

U = a a! . . . o- k _i(P k )o- Q o-i . . . o- k _ 2 {Pk-i) ■ ■ ■ M-PO-Po 

where, for all j < k, Pj is a prefix of o~j(i), for some letter i G A. 

Taking abelianization, since all substitutions have the same incidence ma- 
trix, we get 

l(U) = l(a a 1 ...a k ^i(P k )a o- 1 ...a k ^ 2 (Pk-i)---o-o(Pi)Po) 

= l(a o-i . . . a k -i(Pk)) + 1(<to<ti ■ ■ ■ (7 fc _ 2 (Pfc-i)) + • • • + l(Po) 
M k l(P k ) + M k ' l l{P k ^ + ... Mi(Pi) + MZ(Pi) 

Since the stable and unstable space are invariant under M, the projection 
tt s commutes with M and we get: 

k 

TT s l(U) = J2 Mj Ksl{Pj) 
3=0 

One can always define an adapted norm on the contracting space E s such 
that \MX\ < X\X\, for some A < 1. We obtain : 

M*(tf))l = IEi=o^W(Pj)l 

< Ej=ol^W(Pi)l 

< E"=oA J > s KP,)l 

But there is a finite number of substitutions, hence a finite number of pre- 
fixes, so that \ir s l(Pj)\ is bounded by some constant C; hence the finite sum 
for \7r s l(Pj)\ is bounded, independently of the length of U, by the geometric 
series Zf=o CXj = & 
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This proves that the canonical stepped line associated with the limit point 
stays within bounded distance of the expanding line, and that its projection 
by ir s is bounded. □ 

Definition 3.3. Let M be a positive integral unimodular matrix of the Pisot 
type. Let a be a sequence of substitutions with incidence matrix M. The 
Rauzy fractal associated to the sequence a is the closure of the projection to 
the stable space of the canonical line associated with a limit point of a. It is 
denoted by TZ a . 

We can as above define subsets TZ^ of lZ a corresponding to a particular 
letter; they satisfy a set equation, which relates the Rauzy fractal associated 
with a to the Rauzy fractal associated with Sa. 

Proposition 3.2. We have: 

K= IJ M.TV Sa + nMPt°' j) )) 

Proof. This is an immediate consequence of the equality given above for 
prefixes of the limit point, U = (Tq(U^)Po, by projecting to the stable space 
and taking closure. □ 

We will use this set equation in the next section to prove that Tt a depends 
continuously on the sequence a. 

4. Iterated function systems 

4.1. Hutchinson's theorem. Many fractals are made up of parts that are 
similar to the whole. For example, the middle third Cantor set is the union of 
two similar copies of itself, and the von Koch curve is made up of four similar 
copies. These self-similarities are not only properties of the fractals: they 
may actually be used to define them. Iterated function system do this in a 
unified way and, moreover, often lead to a simple way of finding dimensions. 

A mapping S : W 1 — > M. n is called a contraction if there is a number c 
with < c < 1 such that \S(x) — S(y)\ < c\x — y\ for all x,y € D. Clearly 
any contraction is continuous. If equality holds, i.e. if \S(x) — S(y)\ = 
c\x — y\, then S transforms sets into geometrically similar sets, and we call 
S a contracting similarity. 

Definition 4.1. A finite family of contractions {/i, f'2, ■ ■ ■ , f m }> with m < 2, 
is called an iterated function system or IFS. We call a non-empty compact 
subset K of MP an attractor for the IFS if 

m 

K=(jf i (K). 

i=i 

The fundamental property of an iterated function system is that it deter- 
mines a unique attractor, which is usually a fractal. The proof is a simple 
consequence of Banach fixed point theorem. 
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We denote by Ti the set of nonempty compact subsets of W 1 , and consider 
the map $ defined by: 

$ : U 

771 

K ^(Jfi(K) 

i=i 

Let K a non empty compact set from M n , and x £ W 1 . We define the 
distance d(x, K) = inf yg x d(x, y). We define the Hausdorff distance between 
elements of T~L by 

du{K,K') = max(sup d(x, K), sup d(x,K')). 
Lemma 4.1. [T-L^dy] is a complete space. 

Proof. See [8] □ 
Lemma 4.2. & is a contracting map. 

Proof. If / is a contraction of constant c, then d(f(K),f(K')) < cd(K,K'). 
If {A\, . . . A m ) and (B U ...,B m ) are finite collections of elements of T~L, an 
immediate computation shows that d(L)iAi, Uji?j) < maxj d(Ai, B{). Hence, 
if all the fi are contractions on W 1 of ratio less than c, $ is a contraction on 
% of ratio less than c. □ 

Theorem 4.1. $ admit a unique fixed point. 

Proof. <I> is a contraction of ration c < 1 on a complete space, hence it has a 
unique fixed point by the Banach fixed point theorem. □ 

4.2. Graph-directed Iterated Function Systems. An IFS acts on a 
unique set K, by replacing it by its image [JiLi fi(K), and its fixed point 
is defined by a unique equation. In many cases, the fixed point has a finite 
number of components, and is define by a set of equations, as we saw for 
the Rauzy fractal. This more general framework gives rise to what is called 
Graph-directed Iterated Function Systems (GIFS for short). 

Definition 4.2. Let G be a finite directed graph with set V of vertices 
{l,...,d} and set of edges E. Denote the set of edges leading from i to 
j by Eij. To each e G E associated a contractive mapping f e : K n — > W 1 . If 
for each i there is some outgoing edge we call (G, {/ e }) a GIFS. 

Such a GIFS defines a map $ on H d by 



$ : 




=l..d 

id 



If we define a distance on 7i d by d{(Ki), (K()) = max^d^-, K'^), U d is 
complete for this distance, and $ is a contraction. The same fixed point 
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argument shows that to a GIFS (G,{f e } e eE) there corresponds a unique 
collection of nonempty compacts sets K\ , . . . , K q C W 1 having the property 
that 

K i = U U 

j=l eeEij 

The collection K\, . . . , K q is called GIFS attractor or solution of the GIFS. 

4.3. Fixed point associated to a sequence of IFS. We consider now a 
finite family of k GIFS acting on the same set of vertices V = {1, . . . , d}. To 
each GIFS we associate the contraction &j of ratio Cj < 1 acting on T~L d , as 
defined above, with j G {1, . . . , k}. 

We define C = {1, . . . , A;} N , the set of sequences with value in {1, . . . , k}, 
with the shift acting on C in a natural way. Our goal is to associate to each 
sequence e = eo, ei, . . . G C a family of d compact sets K e = (K Sj i, . . . , K t ^) 
such that K e = & eo Ks e - 

We consider the set FT~L of continuous functions from C to T~L d with the 
uniform distance, defined for any elements F\,F2 G TH by 

d(Fi,F 2 ) = supd n (F 1 (e),F 2 (e)). 

eeC 

Lemma 4.3. FTl is a complete space. 

Proof. Let F n a Cauchy sequences from F%. This mean for all a > , we can 
find N G N such that, for all p,q > N, d(F p , F q ) < a. Hence in particular, 
for any e G C, F n (e) is a Cauchy sequences in %. Since is a complete 
space, it converge to some limit F^e). Hence the sequence F n converges 
pointwise to some function F^. 

Taking a and N as above, we see that, for any n > N, du(F n {e), F OQ (e)) < 
a, hence d(F n , Foo) < a. This proves that the convergence of F n to F^ is in 
fact uniform, hence F^ is continuous, and the space F% is complete. □ 

Lemma 4.4. We consider the map <I> : F% — > F% defined by : = G 

such that G{e) = $ eo (F(Se)). Then & is a contracting map of ratio c < 1. 

Proof. Let Fi and F 2 in TU. 

We know that, for a given e, <5 eo is a contracting map, hence there exist 
c tQ < 1, depending only on eo, such that 

dn^eoF^^^Se)) < cd n (F 1 (Se),F 2 (Se)) 

If we define c = sup i=1 k Ci, we obtain 

supd w ($F 1 (e),$F 2 (e)) < csupd n (F 1 {Se),F 2 (Se)) 

Hence d($(Fi), $(F 2 )) < cd(Fi,F 2 ); since c < 1, we have proved that $ 
is a contraction. □ 

Theorem 4.2. <I> admet a unique fixed point. 
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Proof. Same as above: $ is a strictly contracting map in a complete space 
J-'Ti, hence it admits a unique fixed point. □ 

4.4. Application to generalized Rauzy fractals. We proved in the pre- 
vious section that, if we have a finite set S of substitution with the same 
matrix of the Pisot type, we can associate to any sequence a of substitutions 
in S a generalized Rauzy fractal lZ a , and that this Rauzy fractal satisfies the 
following equation: 

K= U M.K^ + irMPt ' j) )) 

j,k;w k v =i 

We can view (1Z l a ) as a map 1Z : S — > % d , a i— > (Tv^.), where H is the set of 
nonempty compact subsets of the contracting plane of M. Each substitution 
o"o £ S defines a contraction <& CT0 on 7-L d , the contraction ratio being given 
by the contraction ration of the restriction of M to the contracting space. 

We are exactly in the conditions to apply the theorem above, and we 
obtain: 

Corollaire 4.1. Let S be a finite set of substitutions with the same unimod- 
ular Pisot incidence matrix M. To any sequence a € S , we can associate 
a Rauzy fractal lZ a = U d =l lZ % a in the contracting space of M . This Rauzy 
fractal satisfies the relation 

K= U M.K> Sa + *S{Pt h3) )) 

It depends continuously on a for the natural product topology on 5 N . 

5. EXEMPLES 

A classic example is given by the TVibonacci substitution and the flipped 
Tribonacci substitution, i.e., 

(a — > ab ( a — )■ ab 

b — > ac and o~2 \ b — > ca 

c — ^ a I c — > a 



The incidence matrix of o~\ and a% is I 1 I . The dominant eigenvalue 

\0 1 0/ 

satisfies the relation A" 3 — X 2 — X — 1 = 0, hence the name Tribonacci for 
the substitution. 

The Rauzy fractal of the first substitution is a topological disc [1], simply 
connected , while it is a well known fact that the second fractal is not simply 
connected, and has infinitely generated fundamental group, see |14j ; compare 
Figurejl]. We obtain some nice Rauzy fractals for different sequences w n , in 
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Figure 1. The Rauzy fractals of a\ and 02 

particular we can find some examples with a nontrivial but finitely generated 
fundamental group. 




FIGURE 2. The product Rauzy fractal of <j\ and a<i with 
the sequence u = 222211111111111121 



6. ADDITIONAL REMARKS 

6.1. Substitutions with the same matrix. For any a, the set 1Z a con- 
tains a fundamental domain for the lattice T acting on the contracting plane. 
More precisely, the union U 7g r^o- +7 is the contracting plane. The quotient 
TZu/T is a torus, and the map from the limit point to T extends to a contin- 
uous map from the adic system f2 CT to this torus, which is a semi-conjugacy 
between the shift on tt a and a torus rotation. 
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FIGURE 3. The product Rauzy fractals of a\ and a% with 
the sequence u = 2222111111111111212 




FIGURE 4. The product Rauzy fractals of a\ and <j% with 
the sequence u = 1122112211221122 



On would of course like to prove that, if the matrix is unimodular, this 
semi-conjugacy is a measurable conjugacy (in the non-unimodular case, there 
must also be some p-adic component). 
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It seems very likely that the dimension of the boundary of M. a is a measur- 
able function of a which is shift invariant; hence, by ergodicity of the shift 
for any Bernoulli measure, it must be almost everywhere constant for this 
measure. However, this dimension is not constant, since it is generally not 
the same for the various substitutions; one would like to be able to compute 
the generic value for a given Bernoulli measure. 

6.2. The general adic case. It would be much more interesting to ex- 
tend this to the case of substitutions with different matrices, but this poses 
difficult problems. 

We can define the symbolic system for any primitive sequence of substi- 
tutions, but after that we meet several obstructions. 

The first one is to find the equivalent of the Perron-Frobenius eigenvector, 
or the asymptotic direction for the stepped line. This is not always possible: 
the existence of this asymptotic direction is equivalent to the existence of a 
well-defined frequency of the letters, or to the unique ergodicity. But there 
are 5"-adic systems which are not uniquely ergodic, and we can even find 
such systems on 4 letters which are generated by a primitive sequence of 
unimodular substitutions. 

However, such examples seem exceptional, and it might well be possible 
to find quite general conditions which ensure unique ergodicity of the 5-adic 
system. 

But if this preliminary step is satisfied, we still need to find the equivalent 
of the Pisot condition. One should not expect to find an exact equivalent: 
the Pisot condition only makes sense for algebraic numbers, and for most se- 
quences of substitution, the equivalent of the Perron Frobenius eigenvector 
should have transcendant coordinates, and there is no equivalent for eigen- 
values. 

The best condition is probably to ask for a bounded distance between 
the stepped line associated with the limit point and the asymptotic line. 
Another equivalent condition is the condition of bounded balance; we say 
that a sequence is C-balanced if there exists a constant C such that, for any 
two words U, V of same length occuring in the sequence, and any letter a, 
\\U\ a -\V\ a \ < C. 

Any sequence with bounded balance can define a generalized Rauzy frac- 
tal, and it would be very interesting to give sufficiently general condition 
which imply bounded balance; as shown by the case of unbalanced epistur- 
mian sequences, this does not seem to be an easy question. 
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